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Abstract 

In  a  hypersonic  boundary  layer  over  a  wall  of  variable  curvature,  the  region  most  sus¬ 
ceptible  to  Gortler  vortices  is  the  temperature  adjustment  layer  sitting  at  the  edge  of  the 
boundary  layer  (Hall  &  Fu  (1989),  Fu,  Hall  &  Blackaby  (1990)).  This  temperature  ad¬ 
justment  layer  is  also  the  most  dangerous  site  for  Rayleigh  instability  (Cowley  &  Hall 
(1990),  Smith  &  Brown  (1990)  and  Blackaby,  Cowley  and  Hall  (1990)).  In  this  paper,  we 
investigate  how  the  existence  of  large  amplitude  Gortler  vortices  affects  the  growth  rate 
of  Rayleigh  instability.  The  effects  of  wall  cooling  and  gas  dissociation  on  this  instability 
are  also  studied.  We  find  that  all  these  mechanisms  increase  the  growth  rate  of  Rayleigh 
instability  and  are  therefore  destabilizing. 
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1  Introduction 


Curved  surfaces  are  unavoidable  in  many  engineering  designs.  The  engine  inlets  and  the 
control  surfaces  of  hypersonic  vehicles  are  such  examples.  Peculiar  to  curved  surfaces  is  the 
Gortler  instability  mechanism  which  can  be  induced  by  wall  curvature.  In  this  paper,  we 
continue  our  systematic  studies  on  the  Gortler  instability  mechanism  in  the  hypersonic  context. 
The  starting  point  of  our  present  study  is  the  nonlinear  theory  presented  in  our  previous 
paper,  Fu  and  Hall  (1991).  In  that  paper  we  have  shown  that  when  Gortler  vortices  evolve 
downstream  of  the  neutral  position,  a  large  amplitude  vortex  structure  can  be  established 
under  the  combined  action  of  nonlinearity  and  viscosity.  This  structure  consists  of  a  region  of 
vortex  activity  bounded  by  two  viscous  transition  layers  over  which  the  vortices  are  reduced  to 
zero  exponentially.  Above  the  upper  transition  layer  and  below  the  lower  transition  layer  there 
is  only  the  mean  flow.  The  determination  of  the  locations  of  the  two  transition  layers  gives 
rise  to  a  free  boundary  problem  which  we  have  solved  for  different  curvature  distributions. 
After  such  a  large  amplitude  Gortler  vortex  structure  has  been  established,  various  travelling 
waves  may  be  triggered  off  in  the  form  of  secondary  instabilities.  Experimental  results  given 
by  Swearingen  and  Blackwelder  (1987)  suggest  that  at  least  two  forms  of  secondary  instability 
may  exist.  The  first  instability  is  located  at  the  two  transition  layers,  7r/2  out  of  phase  in  the 
spanwise  direction  with  the  existing  Gortler  vortices  and  thus  leads  to  wavy  vortex  boundaries. 
The  other  possible  instability  is  in  phase  with  the  existing  Gortler  vortices  and  has  an  inviscid 
nature.  Thus  it  is  a  Rayleigh  instability.  In  our  previous  paper  Fu  &  Hall  (1991),  the  wavy  type 
of  instability  has  been  luUy  studied  and  we  have  shown  that  a  family  of  neutral  travelling  wave 
modes  may  exist.  In  the  present  paper,  we  shall  investigate  the  Rayleigh  secondary  instability 
which  the  large  Gortler  vortex  structure  may  suffer. 

The  present  secondary  instability  problem  can  also  be  interpreted  in  the  following  way. 
It  has  been  shown  in  Hall  &  Fu  (1989)  and  Fu,  Hall  &  Blackaby  (1990)  that  in  a  hypersonic 
boundary  layer  over  a  wall  of  variable  curvature,  the  region  most  susceptible  to  Gortler  vortices 
is  the  temperature  adjustment  layer  sitting  at  the  edge  of  the  boundary  layer.  This  temperature 
adjustment  layer  is  also  the  most  dangerous  site  for  Rayleigh  instability,  as  has  been  shown 
by  Cowley  &  Hall  (1990),  Smith  &  Brown  (1990),  and  Blackaby,  Cowley  and  Hall  (1990).  It 
is  therefore  of  interest  to  investigate  how  Gortler  vortices  interact  with  Rayleigh  travelling 
waves  in  the  temperature  adjustment  layer.  Of  greater  importance  is  the  question  whether 
the  existence  of  Gortler  vortices  in  the  boundary  layer  makes  the  layer  more  susceptible  to  a 
Rayleigh  instability.  It  has  been  established  in  Fu  &  Hall  (1991)  that  as  Gortler  vortices  develop 
downstream  of  the  neutrally  stable  location,  nonlinear  effects  produce  a  mean  flow  correction 
as  large  as  the  original  basic  state.  Thus  at  locations  0(1)  downstream  of  the  neutral  position, 
the  original  basic  state  is  completely  altered  by  Gortler  vortices.  Our  main  concern  in  this 
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paper  is  to  understand  how  such  alterations  of  the  basic  state  by  Gortler  vortices  aifect  the 
growth  rate  of  Rayleigh  instability. 

The  other  problems  which  we  consider  here  are  the  influences  of  gas  dissociation  and  wall 
cooling  on  Rayleigh  instability  in  the  temperature  adjustment  layer.  The  influence  of  gas 
dissociation  does  not  seem  to  have  been  studied  before  in  the  large  Mach  number  limit,  whilst 
the  influence  of  wall  cooiiug  has  previously  been  examined  by  Blackaby,  Cowley  and  Hall 
(1990),  but  their  analysis  is  for  unit  Prandtl  number.  Here  we  shall  re-examine  this  problem 
for  Prandtl  number  equal  to  0.72  which  is  more  relevant  to  air. 

This  paper  is  organised  as  follows.  In  the  next  section,  we  state  our  problem  and  give 
asymptotic  solutions  for  the  basic  state,  whilst  in  section  3  we  describe  the  large  amplitude 
Gortler  vortex  structure  and  the  resulting  mean  flow.  The  numerical  solution  of  the  latter  is 
then  discussed  in  section  4,  in  preparation  for  the  solution  of  Rayleigh  equation.  In  section 
5,  we  first  formulate  the  Rayleigh  instability  problem  and  then  study  the  influence  of  wall 
cooling  and  gas  dissociation  in  the  absence  of  Gortler  vortices  on  the  Rayleigh  modes.  Finally 
we  investigate  the  influence  of  Gortler  vortices  on  Rayleigh  modes.  We  discuss  our  numerical 
results  and  then  draw  some  conclusions. 


2  The  original  basic  state 


Consider  a  hypersonic  boundary  layer  over  a  rigid  wall  of  variable  curvatura  {Ij A)k{x* j L)., 
where  i  is  a  typical  streamwise  length  scale  and  /I  is  a  lengthscale  characterizing  the  radius  of 
curvature  of  the  wall.  We  choose  a  curvilinear  coordinate  system  (i*,  y*,  z*)  with  x*  measuring 
distance  along  the  wall,  y*  perpendicular  to  the  wall  and  z*  in  the  spanwise  direction.  The 
corresponding  velocity  components  are  denoted  by  (u*,v*,tn*)  and  density,  temperature  and 
viscosity  by  p*,T*'  and  /z*  respectively.  The  free  stream  values  of  these  quantities  will  be 
signified  by  a  subscript  oo.  We  define  a  curvature  parameter  6  by 

6  =  (2.1) 


and  consider  the  limit  6 


0  with  the  Reynolds  number  R  defined  by 


<oLpl 


taken  to  be  large  so  that  the  Gortler  number 


(2.2) 


G  = 


(2.3) 


is  0(1).  In  the  following  analysis,  coordinates  {x*,y*,z*)  are  scaled  on  {L,R~^^^L,R~^/'^L), 
the  velocity  {u*,v*,w*)  is  scaled  on  and  other  quantities  such  as 

p*,T*,  and  /z*  are  scaled  on  their  free  stream  values  with  the  only  exception  that  the  pressure 
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p*  is  scaled  on  /»to“oo  viscosity  is  scaled  on  (Such  a  scaling  is  only  appropriate 

to  the  Gdrtler  problem  and  in  the  Rayleigh  problem  to  be  discussed  in  section  5  a  different 
scaling  will  be  used).  All  dimensionless  quantities  will  be  denoted  by  the  same  letters  without 
a  superscript  *.  Then  the  Navier-Stokes  equations  are  given  by 


Du  dp  d ,  du.  d  .  du 


dy^^dy' '  dz^’^dz^ 


,  1/-.  2  D^P  ,  ^ 

p— — I-  —Gku  — 

^  Dt  2  dy  dy 


^  dxp^^  dy^ 


d  ,  dv.  d  ,  dv. 


iJUJ  Tsdp  d  1  -  .  4.  « 


d  ,  dv0.  d ,  dw.  d ,  dw. 


= a(7  -  + (|)  V  (7  -  -  P(|)7i^ 

d  ,jT.i  d  ,,dT. 
a  dy^^  dy^  a  dz^^  dz^' 

7Af^P  =  (1  +  a)pT.  (2.9) 

Here  we  have  used  a  mixed  notation  in  which  (v].,V2,V3)  is  identified  with  (u,v,w)  and 
(xi,®2,a:3)  with  (x,y,z).  Repeated  suffices  /?  signify  summation  from  1  to  3.  The  functions 
A,  k,  Cp  and  h  denote  in  turn  the  bulk  viscosity,  the  coefficient  of  heat  conduction,  the  specific 
heat  at  constant  pressure  and  the  enthalpy  per  unit  mass.  The  constants  7,  M  and  a  are  in 
turn  the  ratio  of  specific  heats,  the  Mach  number  and  the  Prandtl  number  defined  by 


.V  -  £e«  m2  -  Jk-  -  ^ 


PopCpoo 


where  K  is  the  gas  constant  and  a^o  =  is  the  sound  speed  in  the  free  stream.  Finally, 

the  function  a  in  the  equation  of  state  (2.9)  denotes  the  percentage  by  mass  of  the  mixture 
which  has  been  dissociated.  (We  assume  that  the  gas  is  an  ideal  diatomic  gas,  say  A2.  After 
dissociation  has  taken  place,  each  A2  molecule  dissociates  into  two  A  atoms  and  the  gas 
becomes  a  gas  mixture  of  A2  and  A),  In  equations  (2.5)-(2.8),  the  operator  D/Dt  is  the 
material  derivative  and  it  has  the  usual  expression  appropriate  to  a  rectangular  coordinate 
system. 

The  basic  state  is  given  by 


(u,  V,  w)  =  {u{x,y),  v{x,y),  0),  T  =  f{x,y), 
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P  =  p{x,  y),  P  =  ii{x,  y).  (2.10) 

By  substituting  (2.10)  into  the  governing  equations  (2.4)-(2.9)  it  is  straightforward  to  obtain 
the  reduced  equations  satisfied  by  the  basic  state.  The  reader  is  referred  to  the  book  by 
Stewartson  (1964)  for  a  detailed  discussion  of  these  equations.  If  we  define  the  Howarth- 
Dorodnitsyn  variable  y  and  a  similarity  variable  r)  by 

f  =  /  and  tj  =  -|=,  (2.11) 

Jo  v2* 

then  the  continuity  equation  is  satisfied  if  u  and  v  are  written  as 

u  =  /'(?;),  V  =  r  -dr]].  (2.12) 

V2a;  p  Jo  p 

Here  the  functions  /(q)  and  f{r])  must  satisfy 


//"  +  (W7  =  0,  (2.13) 

kpkf'Y  +  cjf  +  =  0,  (2.14) 

o 

if  the  *-momentum  and  energy  equations  are  to  be  satisfied.  These  equations  must  then  be 
solved  such  that  /,  /'  vanish  at  the  wall,  /',  T  =  1  at  infinity  and  either  T'  =  0  or  T  specified 
at  the  wall.  The  y-moraentum  equation  gives 


to  leading  order  so  that  p  =  p(r).  In  the  following  analysis,  we  assume  that  there  is  no  pressure 
gradient  along  the  streamwise  direction  and  therefore  we  can  takep  =  constant.  Equation  (2.9) 
then  gives 


[l  +  a(f)J^=l, 


(2.15) 


In  the  Gorlter  problem,  the  pressure  perturbation  is  0(l/ii)  relative  to  the  unperturbed  pres¬ 
sure  and  therefore  the  above  relation  for  the  basic  state  can  be  extended  to 


[i-fa(r)]pr=i, 


(2.16) 


which  is  valid  for  the  total  flow. 

We  now  give  a  summary  of  the  asymptotic  solutions  of  the  basic  state  equations  (2.13) 
and  (2,14)  for  an  ideal  dissociating  gas.  A  detailed  derivation  of  these  results  can  be  found  in 
Fu,  Hall  &  Blackaby  (1990).  The  corresponding  results  for  an  ideal  undissocialed  gos  can  be 
obtained  by  putting  a  =  0  in  the  following  equations. 

An  ideal  dissociating  gas  is  a  diatomic  gas  which  satisfies  the  following  dissociation  law: 

v2  « ,  y 


a* 


Pd 


1 


-Ta/T 


a* 


P  Td 


(2.17) 


A 


(see  Becker  (1968),  p.36,  or  Lighthill  (1957),  p.6).  Here  pd  and  Td  are  respectively  the  char¬ 
acteristic  pressure  and  temperature  for  dissociation.  On  rewriting  this  relation  in  terms  of 
nondimensional  vzuriables  and  using  the  fact  that  p  is  a  constant,  we  have 

f%.  Tj 

(2.18) 


,  2  =  £±fe-^, 

1  Poo 


where  pd  =  pdH^d)  is  the  characteristic  density  for  dissociation.  For  the  purpose  of  asymp¬ 
totic  analysis,  it  is  convenient  to  define  two  new  constants  a  and  b  by 


o  =  b  = 


Td 


Poo 


Then  equation  (2.18)  becomes 


a* 


1 


aT 


(2.19) 


(2.20) 


which  displays  the  physical  fact  that  dissociation  will  take  place  in  the  hottest  region  where 

f  =  0(M2). 

We  assume  that  after  dissociation  each  component  of  the  gas  mixture  has  its  viscosity  given 
by  Sutherland’s  law.  We  further  assume  that  the  viscosity  of  the  gas  mixture  is  described  by 
Wilke’s  law  and  the  coefficient  of  heat  conduction  given  by  Wassiljewa’s  formula.  It  can  be 
shown  that  in  the  large  Mach  number  limit,  the  boundary  layer  splits  into  two  sublayers:  a 
wall  layer  and  a  temperature  adjustment  layer.  In  the  wall  layer,  f  M^,  f 

so  it  is  appropriate  to  define  0(1)  variables  Y,t  and  F{Y)  by 


Y  =  M^I\  f  =  M-^f,  F{Y)  =  M^I'^f. 

Then  equations  (2.13)  and  (2.14)  give 


1  4-  m  /  /i2(«)  T' 


M{o^)  {F") 


|//\2 


a 


where 


^  def.,  -^3(1  -  Q)  _ ^2'd3a 

^  2a  +  i43(l-Q)  i43a -I- ^2(1  -  a)’ 


,  def.  ^3(1  -  «) 

ft2  = 


+  “f 


(10^2^3/7)0 


2a  +  i3(l-o)  ^430  +  ^2(1-0)’ 

=  l  +  ~o  +  i(l  + |)". -0(1-02). 


(2.21) 

(2.22) 

(2.23) 

(2.24) 

(2.25) 

(2.26) 


Here  A2,  are  constants  and  m  is  the  constant  appearing  in  Sutherland  viscosity  law  for  the 
undissociated  gas.  On  rewritting  (2.20)  in  terms  of  T,  we  have 

v2 


O* 


1  — 


(2.27) 
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Equations  (2.22)-(2.27)  are  to  be  solved  numerically  subject  to  the  conditions 


F(0)  =  1^(0)  =  0,  r(oo)  =  0,  1^(00)  =  1, 
7^(0)  =  0  if  the  wall  is  thermally  insulated, 


(2.28) 


f'(O)  =  nTfu  if  the  wall  is  under  cooling, 

where  Ty,  is  the  wall  temperature  scaled  on  M^Tco  when  the  wall  is  thermally  insulated  and  n 
is  the  wall  cooling  coefficient. 

As  y  00,  t  0  and  a  decays  to  zero  exponentially.  Prom  (2.24)-(2.26)  we  have 
hi  — ►  1,  ^2  1  and  Cp  — ►  1.  Equations  (2.22)  and  (2.23)  can  then  be  approximated  by 


pw 

(1  +  m)(4=y  +  PF"  =  0, 

VT 

+  Ff-  +  (7  -  1)(1  +  =  0, 


which  have  asymptotic  solution 


s:Y  ~  d  +  — — —  +  •  •  •  r  =  _ _ ! _ + 

^  ’  <r  (y-/3)< 


(2.29) 


(2.30) 


(2.31) 


'  -  ‘  ■  (F-/?)3/.  ■  ’  -  -  L  (r  J  {Y-py  "  ’ 

where  both  D  and  0  are  constants.  In  the  temperature  adjustment  layer  t}  =  0(1),  there  is  no 
dissociation  and  in  order  to  match  with  (2.31),  the  solutions  for  /  and  f  must  expand  as 


Here  Mi  is  defined  by 


^  ^  Ml  ’ 

f  =  t(,/)  +  ..-. 


Mi  = 


(2.32) 

(2.33) 


On  substituting  (2.32)  and  (2.33)  into  (2.13)  and  (2.14),  we  obtain  to  leading  order 


(1  +  m) 


f  +  m 


f"  +  nf"  =  0, 


+  r]f'  =  0. 

O’  \T  +  m  j 

These  two  equations  are  to  be  solved  numerically  subject  to  the  matching  conditions 


as  /(»?)-  ^,  T 


9(1  +  my  1 


and,  at  infinity 


/'(oo)  =  0,  ^(oo)  =  1.. 


(2.34) 


(2.35) 


(2.36) 


(2.37) 
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It  can  be  shown  from  (2.34)  and  (2.35)  that 


/"(»?)  = 


D 


(l  +  m)i+V<^  (12) 


l/<r 


1-1/0 


(2.38) 


so  that  after  (2.35)  has  been  solved  numerically,  the  function  /"(»/)  can  be  computed  easily 
from  this  equation.  Also,  we  note  that  whilst  the  solution  of  (2.35)  is  independent  of  the 
wall  layer  solutions  and  thus  of  the  conditions  at  the  wall,  the  function  /  is  dependent  on 
the  wall  layer  solutions  through  the  matching  constsmt  D.  As  will  be  shown  later,  D  is  an 
important  constant  and  it  is  through  this  constant  that  the  influence  of  wall  cooling  and  gas 
dissociation  affect  the  growth  rate  of  Rayleigh  instability  in  the  temperature  adjustment  layer. 
In  table  1,  we  have  shown  the  dependence  of  D  on  wall  cooling  and  gas  dissociation.  In 
our  numerical  integration  of  the  boundary  layer  equations  (2.22^  and  (2.23),  we  have  taken 
o  =  0.72,  m  =  0.508,  7  =  1.4,  Aj  =  A3  =  1,  o  =  1.225  x  10®,  6  =  3.265.  As  a  check,  we 
have  also  used  our  program  to  obtain  the  value  of  D  when  (7=1,  wall  is  insulated  and  gas 
dissociation  neglected,  for  which  the  exact  solution  for  D  is  possible  and  is  given  by 

7-1 

Our  numerical  solution  yields  D  =  17.02,  whilst  (2.39)  gives  D  =  17.08. 


(2.39) 


Table  1  Dependence  of  D  on  gas  dissociation  and  wall  cooling 


0 

II 

n=0.4 

n=0.6 

II 

0 

bo 

n=l 

Ideal  gas 

395 

313 

257 

216 

186 

Real  gas 

405 

326 

271 

232 

206 

3  The  large  Gortler  vortex  structure 

Because  of  the  curvature  of  the  <vall,  the  hypersonic  boundary  layer  described  in  the  previous 
section  may  lose  stability  to  Gortler  vortices.  In  the  linear  Gortler  instability  analysis,  we  are 
concerned  with  the  determination  of  the  conditions  under  which  Gortler  vortices  grow  in  the 
streamwise  direction.  To  this  end,  we  superimpose  a  steady  periodic  stationary  vortex  structure 
with  wavenumber  a  in  the  spanwise  direction  on  the  basic  state  (2.10)  and  the  perturbation 
equations  are  found  by  linearizing  the  Navier-Stokes  equations  about  the  basic  state.  These 
linear  equations  have  been  fully  discussed  in  our  previous  paper  Fu,  Hall  and  Hlarkaby  (1990).: 
It  was  shown  there  that  for  the  wall  mode  which  has  wavelength  comparable  \m  '  li  the  boundary 
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layer  thickness,  the  neutral  Gdrtler  number  is  a  decreasing  function  of  the  local  wavenumber. 
As  the  latter  increases,  the  centre  of  vortex  activity  moves  towards  the  temperature  adjustment 
layer  and  the  Gdrtler  number  tends  to  a  constant  which  is  the  leading  order  term  in  the  Gdrtler 
number  expansion  for  the  mode  trapped  in  the  temperature  adjustment  layer.  It  is  this  mode 
that  is  most  susceptible  to  Gdrtler  vortices  since  it  has  a  smaller  Gdrtler  number  than  any 
other  mode. 

As  is  typical  of  Gdrtler  vortices  in  growing  boundary  layers,  the  evolution  of  Gdrtler  vortices 
in  the  temperature  adjustment  layer  is  dominated  by  nonparallel  effects.  It  was  shown  in  Fu, 
Hall  and  Blackaby  (1990)  that  in  the  hypersonic  limit  such  nonparallel  effects  operate  mainly 
through  the  curvature  of  the  basic  state  in  the  general  case  when  the  wall  curvature 

is  not  proportional  to  l/(2a;)^/^.  Thus  only  when  the  wavenumber  o  is  as  large  as  of  order 
do  nonparallel  effects  become  negligible  and  the  following  asymptotic  expression  for  the 
neutral  Gdrtler  number  can  be  obtained: 


G  = 


/«(a:„)(2a:„)3/2 


+  gQa^  +  a® 


3ffo  , 

>2  ' 


V^V  2^0 


(3.1) 


where 

B  ‘‘If-  Um  rf{,)d,,  Jo  =  (3.2) 

M~*oo  Jo  aK{Xn)Ti 

and  Tq  =  fi  =  fio  =  jiifo).  The  constant  tj*  denotes  the  centre  of  vortex 

activity  and  has  the  numerical  value  of  3.001  when  a  =  0.72,  m  =  0.509,  whilst  Xn  is  the 
neutral  position.  In  (3.1)  the  first  term  is  due  to  the  curvature  of  the  basic  state  and  other 
terms  are  due  to  viscous  effects.  It  is  clear  that  a  =  is  the  order  at  which  viscous 

effects  become  comparable  with  the  effects  of  centrifugal  acceleration  due  to  the  curvature  of 
the  basic  state. 

Equation  (3.1)  is  a  relation  between  the  ner+ral  Gdrtler  number  G,  the  wavenumber  a  and 
the  neutral  position  a;„.  In  theory  this  relation  can  be  inverted  to  give  an  expression  for  Xn  as 
a  function  of  G  and  a.  Thus  for  a  given  Gdrtler  number  and  a  wavenumber,  we  know  where 
the  vortices  will  become  neutrally  stable.  After  the  neutral  position  has  been  determined,  our 
next  task  is  to  investigate  how  Gdrtler  vortices  will  grow  beyond  the  neutral  position.  This 
task  was  taken  up  in  the  nonlinear  theory  presented  in  our  previous  paper  Fu  and  Hall  (1991). 
It  is  shown  there  that  initially  at  the  neutral  position,  Gdrtler  vortices  are  trapped  in  a  thin 
viscous  layer  of  0(€*/®)  thickness  (due  to  their  small  wavelength  nature);  but  as  they  evolve 
downstream  of  the  neutral  position,  they  spread  into  the  temperature  adjustment  layer  and 
develop  into  a  large  amplitude  vortex  structure.  This  structure  consists  of  a  core  region  of 
vortex  activity  bounded  by  two  viscous  transition  layers  where  Gdrtler  vortices  are  forced  to 
decay  to  zero  exponentially.  The  total  flow  is  written  as 

u  =  u-\-  -77-17,  v  =  v-\-V,  w  =  W^ 

Ml 
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p  —  p+-^(P  +  P)i  T-f  +  T.  (3.3) 

Here  {ii,v,p  +  p/R,f)  is  the  mean  flow;  whilst  {U/Mi,V,W,P,T)  is  the  hannonic  pMt  of 
Gortler  vortices.  The  mean  velocity  components  u  and  v  have  the  following  expressions: 


dfix,t]) 


Here,  in  analogy  with  (2.32)  and  (2.33),  f(x,r})  and  f{x,T))  expand  as 


f{x,v)  V  + 


T  =  Tix,r))  + 


Equation  (3.4)  are  similar  to  (2.12)  except  that  /  and  f  in  (3.4)  are  also  functions  of  x.  This 
is  due  to  the  fact  that  after  the  basic  state  has  been  reinforced  by  nonlinear  interaction,  it 
can  no  longer  be  described  by  the  similarity  variable  r}  along.  Naturally,  if  vg  is  not  present 
in  (3.4b),  we  do  not  expect  u  and  v  given  by  (3.4)  to  satify  the  continuity  equation.  Thus 
the  function  vs(x,  tj)  in  (3.4b)  is  added  so  that  the  continuity  equation  could  be  satisfied.  The 
function  1(f)  in  (3.4b)  is  the  integral  of  the  mean  temperature  from  0  to  77  (thus  the  wall  layer 
temperature  is  also  involved).  The  explicit  expressions  for  it  and  vg  can  be  found  in  Fu  &  Hall 
(1991).  We  do  not  write  them  out  here  since  they  are  not  needed  in  our  following  analysis. 
The  similarity  variable  i]  —  i}(x,  y)  is  defined  by 


1  n 

V2iJo 


dy 

T(x,ri(x,y)y 


where  the  function  f  in  the  integrand  is  understood  to  be  the  composite  solution  of  the  mean 
temperature  (i.e.  the  wall  layer  temperature  and  the  mean  temperature  in  the  temperature 
adjustment  layer). 

In  the  core  region  of  vortex  activity,  the  mean  flow  part  and  the  harmonic  part  of  the  flow 
have  the  following  form  of  solutions: 

=  fo(x,v)  +  (fi{x,v)  + 

f(x,  ri)  =  fo{x,  i)  +  cfi(x,  »?)  +  •••, 
vg  -  Vg(x,r])  +  €Vs{x,r))+---, 


U  =  e{E{U^  +  €Ul  +  -- 

•)  +  €E\UiT- 

•)  +  •• 

•  +  C.C.}, 

y  =  {^(^0^  +  eVi^  +  • 

■ .)  +  ,E\Vo^  +  • 

••)  +  • 

■•  +  c.c] 

W  =  {e(WJ  +  eWl  +  -  • 

■)  +  eE^W^  +  - 

••)  +  ■•' 

•  •  +  c.c] 

P  =  {e(P^  +  ePl  +  . 

...)  +  eE\P^  +  - 

•••)  +  • 

••  +  (7.(7.} 
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T  =  c  +  iB\  +  •••)  +  c£’*(do  +  "•)  +  •••  +  C.C.|’ , 

where 

E  =  exp(— ),  €=-  (3.8) 

€  a 

and  C.C.  denotes  the  conjugate.  It  has  been  shown  in  Fu  and  Hall  (1991)  that  the  mean  flow 
quantities  To,  fo  and  the  fundamental  Vq  satisfies  the  following  three  equations: 


<T  dTo  V2® 

+  Hix)  ~  "• 

(3.9) 

fihdfo'' 

[fodt}) 

1  +  ”  a,  '  9*  -  (x)'“  a,  ’ 

(3.10) 

1  a  ( fioafc]  aT„  afo  2vS 
<r  a,  (5b  a,  J  +  ^  a,  ^  ax  -  a(x) 

(3.11) 

Here  II(x)  is  definec  by 

^^{/c(®n)(2a:„)3/2  (2x)3/2}'^2''^“’ 

(3.12) 

and  fio  =  /i(To).  The  constant  N  in  (3.12)  is  given  by  and  is  of  0(1)  since  we 

have  assumed  that  the  wavenumber  is  of  order  In  our  numerical  calculations  N  is  taken 
to  be  unity. 

Equations  (3.10)  and  (3.11)  are  the  “modified”  forms  of  the  basic  state  equations  (2.34) 
and  (2.35).  The  appearance  of  the  fundamental  in  the  forcing  terms  on  the  right  hand 
sides  reflects  the  fact  that  the  basic  state  is  now  completely  altered  by  nonlinear  interaction. 

The  mean  flow  temperature  To  can  be  first  solved  from  (3.9).  We  note  that  as  x  tends  to 
the  neutral  position  To  tends  to  the  basic  state  temperature  and  (3.9)  then  reduces  to  the 
condition  for  the  Gortler  vortices  to  be  neutrally  stable  at  the  location  t}=  rj*.  The  fact  that 
equation  (3.9)  is  now  an  equation  for  determining  the  mean  temperature  To{x,  i])  shows  that 
the  mean  flow  has  to  adjust  itself  so  that  it  is  neutrally  stable  to  vortices  everywhere  in  the 
temperature  adjustment  layer. 

After  To  has  been  found,  equation  (3.11)  can  then  be  solved  to  determine  the  fundamental 
Vq  (and  hence  Uq  etc.  since  they  are  related  to  Vq).  It  is  easy  to  show  that  the  result  of 
solving  (3.11)  is 


2V^lioToJ  (adri'-Ta  dtl'  ^  dri  dx 


) 


dr]. 


(3.13) 


The  two  locations  where  Vq  vanishes  define  the  centres  7]i{x)  and  rj^ix)  of  the  two  transition 
layers.  Finally,  the  solutions  in  the  core  region  of  vortex  activity  are  complete  if  equation 
(3.10)  can  be  solved  to  give  an  expression  for  fo.  This  has  to  be  done  numerically  and  will  be 
discussed  in  the  next  section. 
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In  the  two  transition  layers  centred  at »/  =  171(2:),  »72(®)>  the  Gdrtler  vortices  are  reduced  to 
zero  expontiaUy  and  so  above  the  upper  transition  layer  and  below  the  lower  transition  layer, 
there  is  only  the  mean  flow.  The  latter  expands  as 

^  =  /o(2:,77)  +  0(c),  f{x,r))  =  fo{x,r})  +  0{e),  (3.14) 

where  /o(a:,  t])  and  Tb(x,  77)  satisfy 


dr]  \ 

%dfQ\ 

,fodr]) 

(3.15) 

ii. 

( /to  3  To 

\ ,  aft 

0  ^0  n 

(3.16) 

crdi] 

Vto  dv  ^ 

-  2x-^  =  0. 
ox 

We  note  that  the  governing  equation  (3.16)  for  fo(x,r))  is  decoupled  from  that  for  /o.  We  can 
therefore  solve  (3.16)  on  [0,»7i],  [772,00)  first  subject  to  the  appropriate  boundary  conditions 
at  77  =  0,  00  and  matching  conditions  at  77  =  771,  772  (in  order  to  match  with  Tq  given  by 
(3.9)).  This  gives  rise  to  a  free  boundary  problem  which  has  been  solved  numerically  in  Fu 
and  Hall  (1991)  for  several  curvature  cases.  The  reader  is  refered  to  that  paper  for  a  detailed 
explanation  of  the  numerical  procedure  which  is  used  to  find  the  locations  of  771,  772  and  the 
mean  temperature  field  in  the  whole  region  0  <  77  <  00. 

In  the  following  section,  we  shall  discuss  the  numerical  solutions  of  (3.10)  and  (3.15)  since 
its  solution  for  /o  and  /q  was  not  given  in  our  previous  paper  Fu  and  Hail  (1991)  and  wiU  be 
needed  in  the  solution  of  Rayleigh  equation  in  latter  sections. 


4  Numerical  solution  for  the  mean  streamwise  velocity 


In  this  section,  we  shall  outline  a  numerical  scheme  which  we  have  used  to  integrate  the  partial 
differential  equations  (3.10)  and  (3.15).  For  convenience,  we  shall  drop  the  subsripts  ‘o  ’  on 
/o,  7b  and  /io  in  these  two  equations.  Thus  our  problem  is  to  integrate 


on  [771,772]  and 


dv  \fdril  'dn 


2x—  =  0 


(4.1) 


(4.2) 


on  (0,771]  and  [772,00).  Here  the  mean  temperature  f  and  f  ,  771  and  772  are  already  known 
from  the  solution  of  the  free  boundary  problem  discussed  in  Fu  and  Hall  (1991). 

The  appropriate  boundary  conditions  are 

W 


/-+0as  77-+ 00,  f  - 


+  •  as  77  0, 


(4.3) 
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and  matching  conditions  are 


fix,T]i)  =  f{x,T)i),  fix,r)2)  =  fix,r)2), 

The  interval  (0,  oo)  is  divided  into  three  sub-intervals: 

Ti  :  (0,  %),  r2  :  Ta  :  (%,  oo). 


Our  aim  now  is  to  integrate  (4.1)  in  the  interval  Fj  and  (4.2)  in  the  intervals  Fi  and  F3  using  a 
marching  procedure  in  the  streamwise  direction.  At  each  x,  we  iterate  on  the  values  of  /(x,  7/1) 
and  f{x,r}2)  so  that  the  matching  conditions  can  be  satisfied. 

For  the  purpose  of  our  numerical  calculations,  it  is  necessary  to  work  with  fixed  boundaries 
so  in  Fi ,  F2  and  F3  we  make  the  transformations 

T7  =  7/i(x)e^,  v  =  m  +  Cim-vi)y  n  =  (4.5) 


respectively,  so  that  the  intervals  Fi ,  F2  and  F3  now  become 

F?:(-oo,0),  F<:(0,1),  F^(l,oo). 


(4.6) 


The  additional  exponential  stretching  in  (4.5)  is  introduced  to  accomodate  the  rapid  change 
of  /  near  =  0  (as  indicated  by  (4.3)). 

In  terms  of  the  new  variables  C  &nd  the  governing  equations  for  the  three  intervals 
become 

.«/  .....  ;.«V  _  „  .  i 


Here 


dx 


2*^  +  A2(x,<^,/)^  =  H2(x,C,/,^)i^  in  F^, 


=. 

dCd(: 


.  ,  Vt  e-2^ 

Ai  =  -(1  +  m) 


T  +  m  Vi 

^  g-20 

(1  -f-  m)  •  - 

Vi  T  +  m  m 


+(1  +  m) 


.-20 


m-f 


•(^), 

^1  2\/r(t  -I-  m)2  ’ 


(4.7) 

(4.8) 

(4.9) 


2x 


V2-m 


Wi  +  Civ2-v\)) 


As  =  -(l  +  m) 


y/f 


1 


T  +  m  vi 


„  /,  .  2xtjL\  ,  ,  V  1  m-f  df 

^~\  Vi)'^  Vl  2\^(t  +  m)2  di}’ 


(4.10) 


Equations  (4.7)-(4.9)  are  parabolic  partial  differential  equations,  so  their  solutions  can  be 
obtained  by  a  marching  procedure.  We  shall  now  use  the  solution  of  (4.7)  as  an  illustrative 
example  to  explain  our  numerical  scheme.  K  the  value  of  /  is  known  at  a;  =  x,  then  the 
following  scheme  is  used  to  determine  this  function  at  x  +  c: 


2x  • 


+  i4i(x,<^,-,/,) 


-2/r+/r-i 

/i2 


T)  t  a  A  t  ^ \  /t+1  fi—\ 


2h 


(4.11) 


where  h  is  the  vertical  grid  spacing,  a  superscript  ‘n’  signifies  evaluation  at  new  the  position 
X  +  €  and  a  subscript  signifies  evaluation  at  the  indicated  vertical  grid  point.  If  we  replace  -oo 
by  (^0  and  use  n  mesh  points  in  the  rj  direction,  we  have 


=  <^  4-  ih,  (f>n  =  <f>Q  +  nh  =  0. 


Application  of  (4.11)  for  t  =  1, 2,  •  •  • ,  n  -  1  gives  a  triagonal  system  of  equations  which  can  be 
solved  after  the  following  boundary  conditions  are  applied: 

ZD 


/o  =- 


a{r)ie^y+^/< 


(4.12) 


where  the  value  of  F\  is  a  guess.  Then  the  left  derivative  of  /  at  77  =  j/i  can  be  calculated  by 
using  the  formula 


U.\  -  i!Zl  -  /n”t|-/n"-l 

577  ®  2/1771 


(4.13) 


Similarly,  by  making  a  guess  for  /  (or  f)  aX  t}  =  772,  say  F2,  we  can  solve  (4.8)  and  (4.9)  and 
calculate  the  “right”  derivative  of  /  at  77  =  771,  the  “left”  derivative  of  /  at  77  =  772  and  the 
“right”  derivative  of  /  at  77  =  772.  By  defining  two  functions  Gi  and  G2: 

Gi(Fi,f2)  =  |^U.-|^'u„ 

G2{Fi,F2}  = 


13 


we  can  iterate  on  JFi  and  F2  with  the  aid  of  the  two  dimensional  version  of  the  Newton- Raphson 
method  until  G\  and  G2  become  sufficiently  small  (ensuring  that  the  matching  conditions  (4.4) 
are  satisfied). 

The  above  procedure  shows  how  to  march  the  values  of  /  and  /  one  step  forward  along  the 
streamwise  direction  at  a  given  downstream  location.  The  scheme  is  complete  if  the  values  of 
/  and  /  are  known  at  a  certain  initial  position  x  =  xq.  Su-h  values  are  provided  by  the  weakly 
nonlinear  theory,  as  we  show  below. 

For  small  x  —  Xn,  Gortler  vortices  are  described  by  the  weakly  nonlinear  theory  given  in 
our  previous  paper  Fu  and  Hall  (1991).  It  is  found  there  that  the  mean  streamwise  velocity  in 
the  core  region  of  vortex  activity  may  be  written  as 


(4.14) 


where  /  is  the  same  function  as  that  appearing  in  (2.32),  whilst  Umo  satisfies  the  evolution 
equation 

(2x„fo  d  an  _  2f''  d{Vof 
\  /io  dx  d<i>^)  d(t> 

where  the  scaled  variables  ((^,i)  are  related  to  {x,r})  by 

X-Xn 


(4.15) 


A 

V  —  /I  * 


X  = 


(4.16) 


fl/l  e 

and  To  =  /io  =  /^(Tb)  fd  =  /"{t}*).  In  the  limit  x  -+  00,  it  is  established  in  Fu  and  Hall 

(1991)  that  (Fo)^  has  the  following  asymptotic  solution: 


{Vo?  =  0  •  •  (1  +  2fc)  •  X .  (C^  -  e^)  +  0(x°), 


(4.17) 


(4.18) 


where 

.  m  <f> 

V 

and  a,b,k  are  constants.  The  function  (Fo)*  vanishes  at  f  =  ±C  which,  through  (4.16)  and 
(4.18),  gives  the  locations  >/  =  »/!,%  of  the  two  transition  layers  in  the  small  x  -  x„  limit..  It 
can  then  be  deduced  that  the  corresponding  asymptotic  solution  for  Umo  is  of  the  form 


ttmO  =  +  0(x^/*). 

On  substituting  (4.19)  into  (4.15)  and  equating  the  coefficients  of  we  obtain 

dC^  2  dC  2^ 


where 


H{0  =  H{0,  C  = 


Ixnfob 

2ap.o 


A=-(l  +  2fc) 


6afd 


3/2 


2Xn(T^Td  \2XnTo, 


(4.19) 


(4.20) 


(4.21) 
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It  can  be  shown  that  the  general  solution  of  (4.20)  which  is  bounded  at  ^  =  0  is  given  by 


F(C)  =  AC  +  ci(C"  +  6C), 


(4.22) 


where  the  constant  ci  is  determined  by  matching  (4.22)  with  the  solutions  in  the  two  regions 
above  the  upper  transition  layer  and  below  the  lower  transition  layer.  In  the  region  above  the 
upper  transition  layer,  there  are  no  Gortler  vortices.  The  expansion  (4.19)  is  retained  but  now 
H{Q  satisfies 


d(‘  '^2  dC  2^ 


The  solution  which  is  bounded  as  ^  oo  is 


(4.23) 


j(0  =  c,exp(-^)t;(^,;^),  (4.24) 

where  C2  is  a  disposable  constant  and  U  is  &  parabolic  cylinder  function.  By  matching  (4.22) 
with  (4.24)  at  ^  =  C,  it  is  easy  to  show  that  c\  and  C2  are  given  by 

^  "(3  +  c!)C2+^6(l  +  C3)’  8"^  +  6Cc)}  IV{-,  ^),  (4.25) 


where  and  C3  have  the  expressions 


C3  = 


4  .  _  Ixnfob 

^^"V2d/2o^' 


(4.26) 


We  expect  that  H{0  is  an  odd  function  of  C  and  so  its  solution  for  —00  <  ^  <  -C  can  be 
obtained  from  (4.24)  by  replacing  C  by  and  C2  by  -C2,  respectively.  On  collecting  the  above 
results  together,  we  have 

f  C2exp(-^)tA(|,  ^),  Cc  <  C  <  00 

^(0  =  ^(0  =  I  AC  +  ci(C"  +  6C),  -Cc<C<Co  (4.27) 

[  -C2exp(-^)(^(|.-:^),  -00  <  C  <  Cc 

On  substituting  (4.19)  back  into  (4.14)  and  making  use  of  (4.16),  we  obtain  the  expression 

^  =  +  +  (4.28) 

which  is  valid  for  a;  -  <  1.;  Evaluating  the  first  two  terms  on  the  right  hand  side  of  (4.28)  at 

some  location  near  a;„  then  gives  the  appropriate  initial  conditions  for  the  partial  differential 
equations  (4.7)-(4.9).; 

In  Fig.l,  we  have  shown  the  evolution  of  the  mean  streamwise  velocity  component  down¬ 
stream  of  the  neutral  position  x„  =  0.3.  The  numerical  calculation  corresponds  to  the  curvature 
case  k{x)  =  2x  and  to  an  insulated  wall  without  dissociation.  The  initial  condition  is  obtained 
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by  evaluating  (4.28)  at  xg  =  0.301.  It  is  clear  from  the  graph  that  alteration  of  the  basic 
state  by  nonlinear  Gortler  vortex  interaction  mainly  taJkes  place  in  the  region  of  Gortler  vortex 
activity  and  as  Gortler  vortices  develop  downstream,  they  reinfore  the  basic  state  by  making 
the  mean  streamwise  velocity  grow  monotonically.  In  Fig.2,  we  have  shown  the  effects  of  wall 
cooling  and  gas  dissociation  on  the  development  of  the  mean  streamwise  velocity  by  considering 
three  cases:  insulated  wall  with  dissociation  neglected  (case  I),  insulated  wall  with  dissociation 
taken  into  account  (case  II),  and  cooled  wall  (n=0.6)  with  dissociation  neglected  (case  III). 
The  wall  curvature  and  the  initial  condition  are  the  same  as  in  Fig.l.  Initially  at  the  neutral 
position  xg  =  0.3,  the  amplitude  corresponding  to  case  II  lies  between  those  for  cases  I  and 
III  (since  the  values  of  D  have  this  property  and  the  amplitude  is  proportional  to  D).  Fig.2 
shows  that  this  property  is  stiU  preserved  after  Gortler  vortices  have  evolved  a  distance  of  0.3 
downstream  of  the  neutral  position. 


5  Rayleigh  secondary  instability 


After  all  of  the  mean  flow  quantities  associated  with  the  large  amplitude  Gortler  vortex  struc¬ 
ture  have  been  determined,  we  are  now  in  a  position  to  study  the  Rayleigh  secondary  instability 
mentioned  in  the  introduction.  In  order  to  determine  whether  Rayleigh  secondary  instability 
can  occur  we  superimpose  a  Rayleigh  travelling  wave  structure  on  the  existing  Gortler  vortex 
structure.  Thus  the  total  flow  is  now  written  as 
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Here  the  three  terms  represent  in  turn  the  mean  flow,  the  harmonic  part  of  Gortler  vortices 
and  the  Rayleigh  travelling  wave  structure.  The  small  parameter  6  is  introduced  to  facilitate 
linearization;  whilst  the  short  length  scale  X  and  the  short  time  scale  r  characteristic  of 
Rayleigh  instability  are  defined  by 


X  =  r  = 


(5.2) 


On  substituting  (5.1)  into  the  Navier-Stokes  equations  (2.4)-(2.9)  and  performing  the  usual 
linearization  with  respect  to  S,  we  obtain 


ia(U  -  c)p  -I-  iap*u  + 


d{p*v) 

dy 


d{p*w) 

dz 


=  0, 


(5.3) 
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(5.4) 


ia{U  -  c)u  +  v-^  +  w-^  = - tap, 

ay  02  p* 

.  ..  1  dp 

ta{U-c)v=  — 

p*  ay 

.  .f  .  _  I  dp 

ta{U  -  c)d  +  v~  +  =  - — ~ — toi{U  -  c)p, 

jM^p  =  p*9  +  0*p, 


where 


d<p  dip  M 

U  =  u  +  U,  p*  =  p  +  p,  0*  =  T  +  0.  (5.9) 

We  note  from  (2.16)  that  in  the  temperature  adjustment  layer  where  a  =  0,  p*  and  0*  satisfy 
the  relation  p*0*  =  1.  By  using  this  relation,  we  can  eliminate  u,  v,  w  and  0  in  favour  of  p  in 
(5.3)-(5.8)  and  obtain  the  following  single  equation  for  the  pressure  p: 


'ilxil 

dy^  dz^ 


-  a*  j  p  +  a 


”(P  -  c)’M“  ■  1 


-  — f 

U-c 


dy  dy  "**  dz  dz 


)  ^  0*  \dy  dy  dz  dz  j 


(5.10) 


When  Gortler  vortices  are  absent,  this  equation  reduces  to  the  well-known  Rayleigh  equation 
(see,  for  example,  Cowley  and  Hall  (1990)) 

In  view  of  (5.9),  (3.3)-(3.5)  and  (3.7),  the  functions  U  and  0*  in  (5.10)  have  the  expressions 


0*  =  io(a;,  ?/)+•••  +  (E{0q  -f  •  •  •)  -H - C.C.. 

in  the  region  -^1  <  t)  <  and  have  the  expressions 

0*  =  Tq{x,  ?/)+•••  -I-  c^^^E{0qi  -f  • '  •)  -t-  •  •  •  -f  C.C., 

in  the  regions  -oo  <  rj  <  rji  and  Tf2  <  1]  <  oo. 

We  now  look  for  the  following  form  of  solutions  for  (5.10): 


(5.12) 


(5.13) 


p(x,y,z)  =  Poix,T])  +  €^Pi{x,r],z)  + 
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(6.14) 


By  assuming  that  the  leading  order  z-dependent  term  in  the  expression  for  p  to  be  of  order 
we  will  be  able  to  coniine  our  attention  to  the  discussion  of  Rayleigh  instability  associated  with 
the  generalized  inflexion  point  of  the  mean  flow  and  therefore  exclude  the  possible  instability 
associated  with  the  generalized  inflexion  points  in  the  z-direction  due  to  the  presence  of  Gortler 
vortices. 

On  substituting  (5.12)— (5.14)  into  (5.10)  and  keeping  only  those  terms  which  are  of  order 
we  find  that  Po{x,  rj)  satilies 

in  the  region  —Tji<T]<rj2  and 


d^Po  2  dfodPo  Oj.^2fi2p  _n 
- 5 - ~~a — A - —  0) 

dri^  fo-cdr)  drj 


(5.16) 


in  the  regions  -oo  <  t]  <  rji  and  r)2  <  rj  <  oo.  An  investigation  of  the  disturbance  equations 
in  the  transition  layers  shows  that  these  layers  are  passive  and  that  Pq,  Pqt,  are  therefore 
continuous  at »;  =  »?ij  >12' 

In  the  remaining  part  of  this  section,  we  shall  solve  these  equations  for  different  mean  flows 
to  determine  the  influence  of  the  presence  of  Gortler  vortices,  wall  cooling  and  gas  dissociation 
on  the  growth  rate  of  Rayleigh  modes. 


5.1  The  influence  of  wall  cooling  and  gas  dissociation 

In  the  absence  of  Gortler  vortices,  /o  and  /o  reduces  to  dfo/di]  and  d/o/dr?  reduces  to 
defined  in  (2.38),  whilst  Tq  and  To  reduce  to  f  determined  by  (2.35).  Equation  (5.15) 
then  becomes 

^  -  {V2mfPPo  =  0.  (5.17) 

dr}^  f  -c  or] 

This  equation  has  been  solved  before  by  Blackaby,  Cowley  and  Hall  (1990)  for  unit  Prandtl 
number  without  taking  gas  dissociation  into  account.  Here  we  shall  solve  this  equation  for 
Prandtl  number  a  =  0.72  which  is  more  relevant  to  air.  We  shall  also  study  how  gas  dissociation 
affects  the  growth  rate. 

As  we  have  remarked  at  the  end  of  section  2,  properties  in  the  wall  layer  affect  properties 
in  the  temperature  adjustment  layer  only  through  the  matching  constant  D  which  appears 
in  (2.31).  Thus  gas  dissociation  and  wall  cooling  affect  Rayleigh  instability  simply  through 
changing  the  values  of  D.  We  note  from  (2.38)  that  the  functions  f" ID  and  j' /D  are  inde¬ 
pendent  of  wall  cooling  and  gas  dissociation.  It  is  easily  seen  from  (5.17)  that  if  cq  and  Do 
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denote  the  values  of  c  and  D  when  the  wall  is  insulated  and  gas  dissociation  is  neglected,  then 
for  any  other  situations  the  value  of  c  is  given  by 


c  = 


(5.18) 


In  Table  1,  we  have  given  the  numerical  values  of  D  when  the  wall  is  cooled  or  gas  dissociation 
is  taken  into  account.  It  is  clear  that  both  gas  dissociation  and  wall  cooling  increase  the 
values  of  D  and  thus  through  (5.18)  they  increase  the  growth  rate  of  Rayleigh  instability. 
The  D  values  given  in  Table  1  impUes  that  if  the  wall  is  cooled  to  one  fifth  of  its  insulated 
temperature,  then  the  growth  rate  will  be  more  than  doubled.  Whilst  the  D  values  in  Table  1 
shows  relatively  smaller  infiuence  on  the  growth  rate  by  taking  gas  dissociation  into  account, 
we  should  note,  however,  that  those  values  only  correspond  to  one  set  of  parameters  associated 
with  gas  di'jsociation  (given  in  the  paragraph  above  (2,39)).  In  theory  the  degree  of  gas 
dissociation  can  be  “tuned”  by  choosing  different  values  for  these  parameters  in  our  numerical 
calculation.  In  particular,  the  degree  of  gas  dissociation  can  be  increased  by  choosing  a  larger 
value  for  a  and  and  a  smaller  value  for  h  in  (2.20).  By  (2.19),  this  is  effectively  to  choose  a 
larger  value  for  the  Mach  number  if  the  gas  properties  are  fixed.  Thus  gas  dissociation  may 
have  a  more  significant  effect  on  the  growth  rate  than  that  implied  by  Table  1. 

Equation  (5.17)  has  been  solved  by  using  a  fourth  order  Runge-Kutta  method  to  determine 
the  relation  between  the  growth  rate  (defined  as  V^a  times  the  imaginary  part  of  c)  and  the 
local  wavenumber  y/^a.  The  boundary  conditions  used  are 


Po 


( P  at  dr)),  as  rj  0,  and  Po  ~  exp{-y/2xar]),  as  7/  — *■  oo.  (5.19) 


The  results  from  such  numerical  integration  are  shown  in  Fig.3  and  on  the  same  plot  we 
have  also  shown  Blackaby,  Cowley  and  Hall’s  result  for  a  =  1.  Both  curves  correspond  to  an 
insulated  wall  with  gas  dissociation  neglected.  We  can  see  that  when  a  =  0,72  the  growth 
rate  attains  its  maxium  value  of  2.0735  at  y/txa  =  0.034757.  The  maximum  growth  rate  for 
<T  =  0.72  is  roughly  four  times  as  large  as  the  maximum  growth  rate  corresponding  to  a  =  1.; 
This  is  not  surprising  since  the  growth  rate  is  proportional  to  D  (the  proportionality  factor 
being  a  function  of  <t  as  well  as  the  wavenumber);  the  exact  D  value  for  ct  =  1  given  by  (2.39) 
when  the  wall  is  insulated  is  only  17  when  m  =  0.509,  7  =  1.4,  whilst  the  corresponding  D 
value  for  tr  =  0.72  is  186.  In  Fig.5  and  Fig.6,  we  have  shown  the  real  and  imaginary  parts 
of  the  pressure  for  a  —  0.72  for  a  range  of  wavenumbers.  We  note  that  whilst  the  imaginary 
part  of  the  pressure  decays  to  zero  rapidly  as  rj  —*■  0  ot  t}  -*  00,  the  real  part  of  the  pressure 
decays  more  and  more  slowly  towards  the  edge  of  the  boundary  layer  as  the  local  wavenumber 
k  (=  y/2xa)  decreases.  These  properties  can  easily  be  infered  from  the  decay  relations  (5.19). 

In  order  to  determine  the  neutral  wavenumber,  we  first  note  that  the  generalized  inflexion 
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point  fj  is  located  where 


t 


(5.20) 


Obviously,  it  is  independent  of  wall  cooling  or  gas  dissociation.  A  simple  numerical  calculation 
shows  that  fj  =  1.958.  Then  the  neutrally  stable  Rayleigh  mode  has  c  given  by 


c  =  f'{fj)  =  -18.56. 


(5.21) 


The  corresponding  a  value  is  obtained  by  integrating  (5.17)  rver  0  <  7/  <  oo  with  an  appro¬ 
priate  treatment  at  the  generalized  inflexion  point.  Such  a  calculation  shows  that  the  neutral 
value  of  a  is 

V^a  =  0.286.  (5.22) 

Finally,  in  Fig.4  we  have  shown  the  eflfects  of  gas  dissociation  and  wall  cooling  on  the  growth 
rate  in  the  absence  of  Gortler  vortices.  The  curve  for  an  insulated  wall  with  gas  dissociation 
neglected  is  the  same  as  that  appearing  in  Fig.3,  whilst  the  other  two  curves  are  obtained  by 
using  the  relation  (5.18)  and  Table  1. 


5.2  The  influence  of  Gortler  vortices 

As  Gortler  vortices  develop  downstream  of  the  neutral  position  Xn,  nonlinear  interactions 
produce  an  0(1)  mean  flow  correction  which  completely  alters  the  original  basic  state.  The 
leading  order  mean  streav/ise  velocity  /o,  /o  and  temperature  Tq,  Tq  in  the  Rayleigh  equations 
(5.15)  and  (5.16)  are  calculated  for  any  given  x  by  using  the  numerical  procedures  outlined 
in  section  4  and  our  previous  paper  Fu  and  Hall  (1991),  respectively.  The  Rayleigh  equation 
(5.15)  and  (5.16)  are  then  solved  to  determine  the  dependence  of  the  growth  rate  on  x  (and 
thus  on  Gortler  vortices).  Such  results  for  the  curvature  case  k{x)  =  2x  are  shown  in  Fig.7, 
Fig.8  and  Fig.9,  which  in  turn  correspond  to  an  insulated  wall  with  gas  dissociation  neglected,  a 
cooled  wall  (n=0.6)  with  gas  dissociation  neglected,  and  an  insulated  wall  with  gas  dissociation 
taken  into  account.  In  each  of  these  flgures,  x  =  0.3  is  the  neutral  position  where  there  is  no 
Gortler  vortex  influence.  It  is  clear  that  the  growth  rate  is  increased  by  the  alteration  of  the 
basic  state  by  Gortler  vortices  and  as  Gortler  vortices  develop  downstream,  the  growth  rate 
increases  raonotonically.  We  have  also  considered  walls  with  curvatures  given  by  k(x)  =  (2x)^/* 
and  k(i)  =  For  k(x)  =  {2xf/^,  we  obtained  similar  figures  to  Fig.7-Fig.9,  whilst  for 
k(x)  =  we  found  that  for  small  x  -  the  growth  rate  is  increased  by  Gorter  vortices 
as  the  latter  evolve  downstream,  but  such  increase  is  hardly  noticeable  graphically;  as  Gortler 
vortices  travel  further  downstream,  all  growth  rate  curves  become  coalesced  into  a  single  curve. 
This  is  because  k(x)  =  \/^  is  a  .special  curature  case  in  the  sense  that  at  large  distances 
downstream  of  the  neutral  position,  a  similarity  solution  exists  in  which  the  boundaries  of 
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the  region  of  vortex  activity  and  the  mean  flow  quantities  ail  become  independent  of  x  (and 
depend  only  on  rj,  see  Fu  and  Hall  (1991)). 

5.3  Conclusions 

In  summary,  what  we  have  shown  in  this  paper  is  that  the  presence  of  Gorlter  vortices,  wall 
cooling  and  gas  dissociation  all  have  a  destabilizing  effect  on  Rayleigh  modes.  Although  the 
influence  of  Gorlter  vortices  is  not  so  pronouced  in  a  small  neighbourhood  of  the  neutral 
position  as  that  of  wall  cooling  or  gas  dissociation,  their  existence  in  a  hypersonic  boundary 
layer  can  nevertheless  make  the  later  more  susceptible  to  Rayleigh  modes  if  they  are  allowed 
to  travel  far  enough  distances  from  the  neutrally  stable  position.  The  main  effect  of  Gortler 
vortices  is  to  increase  the  unstable  band  of  Rayleigh  waves,  thus  the  presence  of  large  amplitude 
vortices  is  likely  to  cause  the  boundary  layer  to  become  more  receptive  to  transition  induced 
by  Rayleigh  modes. 
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The  mean  streomwise  velocity  ot  x  =  0.8( 


Fig.  2 


Figure  2.  Profiles  of  the  mean  streamwise  velocity  after  Gortler  vortices  have  evolved  a 
distance  of  0.5  downstream  of  the  neutral  position  Xn  =  0.3  over  a  wall  with  curvature  given 
by  /«(®)  =  2x. 


The  growth  rote 


Fig.3 


Figure  3.  A  comparison  of  the  growth  rate  corresponding  to  (7  =  1  and  0.72  when  Gortler 
vortic^  are  absent. 
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The  real  part  of  the  pressure 
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Figure  5.  The  real  part  of  the  pressure  at  different  values  of  the  wavenumber  {k  =  v'Sxo) 
when  Gortler  vortices  are  absent,  wall  is  insulated  and  gas  dissociation  is  neglected. 
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Figure  6.  The  imaginary  part  of  the  pressure  at  different  values  of  the  wavenumber  {k 
y/^a)  when  Gortler  vortices  are  absent,  wall  is  insulated  and  gas  dissociation  is  neglected. 
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The  growth  rate 


Fig. 8 


Figure  8.  Effects  of  Gortler  vortices  on  the  growth  rate  of  Bayleigh  instability.  The  wall  is 
cooled  (n=0.6)  with  curvature  given  by  k(x)  =  2x  and  gas  dissociation  is  neglected. 
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The  growth  rote 


Fig. 9  Effects  of  Gortler  vortices  (with  dissociation) 


Figure  0.  Effects  of  Gortler  vortices  on  the  growth  rate  of  Rayleigh  instability.  The  wall  is 
insulated  with  curvature  ^ven  by  /c(i)  =  2i  and  gas  dissociation  is  taken  into  account. 


31 


REPORT  DOCUMENTATION  PAGE 

Form  Approved 

0MB  No  0704-0188 

pjDi  -  Durden  ^o'’  th  ^  cc'eaion  ot  information  estimateO  to  average  i  nour  oer  response  induoing  the  time  tor  reviewing  instru«ions  searching  emstmg  data  sources 

gatt'efing  and  maintaining  the  data  needed  and  completing  ana  reviewing  th**  collea-on  of  information  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
coiieaion  of  n’-'ration.  nciudmg  suggestions  for  reducfng  this  ourden  to  Washington  HeaoQuarters  Services  Directorate  for  information  Operations  and  Reports,  1215  Jefferson 
OavisHigpwa,  >jiie ’2C4  Arlington  vA  22202-4302.  and  to  the  Off^e  of  Mar^agement  and  Budget  Paperwork  ReductionProject{0704-0t88)  Washington  DC  20503 

1.  AGENCY  USE  ONLY  (Leave  blank)  2  REPORT  DATE  3.  REPORT  TYPE  AND  OATES  COVERED 

December  1991  Contractor  Report 

4.  TITLE  AND  SUBTITLE 

EFFECTS  OF  GOrTLER  VORTICES,  WALL  COOLING  AND  GAS 
DISSOCIATION  ON  THE  RAYLEIGH  INSTABILITY  IN  A 

HYPERSONIC  BOUNDARY  LAYER 

5.  FUNDING  NUMBERS 

C  NASl-18605 

WU  505-90-52-01 

6.  AUTHOR(S) 

Yibin  Fu 

Philip  Hall 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADORESS(ES) 

Institute  for  Computer  Applications  in  Science 
and  Engineering 

Mail  Stop  132C,  NASA  Langley  Research  Center 

Hampton,  VA  23665-5225 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

ICASE  Report  No.  91-87 

9.  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  ADORESS(ES) 

National  Aeronautics  and  Space  Administration 

Langley  Research  Center 

Hampton,  VA  23665-5225 

10.  SPONSORING /MONITORING 

AGENCY  REPORT  NUMBER 

NASA  CR-189577 

ICASE  Report  No.  91-87 

11.  SUPPLEMENTARY  NOTES 

Langley  Technical  Monitor:  Michael  F.  Card  Submitted  to  Journal  of 

Final  Report  Fluid  Mechanics 

124.  DISTRIBUTION /AVAILABILITY  STATEMENT 

Unclassified  -  Unlimited 

Subject  Category  34  , 

12b  DISTRIBUTION  CODE 

13.  A9STRACT  (MaMimum  200  words) 

'  In  a  hypersonic  b9undary  layer  over  a  wall  of  variable  curvature,  the  region  most 
susceptible  to  Cartier  vortices  is  the  temperature  adjustment  layer  sitting  at  the 
edge  of  the  boundary  layer,  (Hall  &  Fu  (1989),  Fu,  Hall  &  Blackaby  (1990) ) This 
temperature  adjustment  layer  is  also  the  most  dangerous  site  for  Rayleigh  Instabil¬ 
ity.  (Rowley  &  Hall  (1990),  Smith  &  Brown  (1990)  and  Blackaby,  Cowley  and  Hall 
(1990)).  iA  this  paper,  we  investigate^ how  the  existence  of  large  amplitude  Cutler 
vortices  affects  the  growth  rate  of  Rayleigh  instability.  The  effects  of  wall  cool¬ 
ing  and  gas  dissociation  on  this  instability  are  also  studied.  We  find  that  all 
these  mechanisms  increase  the  growth  rate  of  Rayleigh  instability  and  are  therefore 
destabilizing. 
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